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ABSTRACT 



A Near-infrared (1.18 — 1.35 /im) high-resolution spectrum of the 
gravitationally-lensed QSO APM 08279+5255 was obtained with the In- 
■ fraRed Camera and Spectrograph mounted on the Subaru Telescope using the 

^ ■ adaptive optics system. We detected strong Na I D AA 5891, 5897 doublet ab- 

ff^ ■ sorption in high-redshift damped Lya systems (DLAs) at Zabs = 1-062 and 1.181, 

confirming the presence of Na I, which was first reported for the rest-frame UV 



O ■ Na I AA 3303.3, 3303.9 doublet by Petitjean et al. This is the first detection 

• of Na I D absorption in a high-redshift (z > 1) DLA. In addition, we detected 

a new Na I component in the Zabs = 1-062 DLA and four new components 
2 ■ in the Zabs = 1-181 DLA. Using an empirical relationship between Na I and 

■ Hi column density, we found that all components have large H I column 

density (log Nhi [cm~^] > 20.3), so that each component is classified as DLA 



absorption. We also detected strong Na I D absorption associated with a Mg ll 
system at Zabs = 1.173. Because no other metal absorption lines were detected 
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in this system at the velocity of the Na I absorption in previously reported 
optical spectra (observed 3.6 years ago), we interpret this Na I absorption 
cloud probably appeared in the hne of sight toward the QSO after the optical 
observation. This newly found cloud is likely to be a DLA based upon its 
large estimated H I column density. We found that the NNai/Ncaii ratios in 
these DLAs are systematically smaller than those observed in the Galaxy; they 
are more consistent with the ratios seen in the Large Magellanic Cloud. This 
is consistent with dust depletion generally being smaller in lower metallicity 
environments. However, all five clouds of the Zabs = 1.181 system have a high 
NNai/Ncaii ratio, which is characteristic of cold dense gas. We tentatively suggest 
that the host galaxy of this system may be the most significant contributor to 
the gravitational-lens toward APM 08279+5255. 

Subject headings: galaxies: formation-gravitational lensing-intergalactic medium- 
quasars: absorption lines-quasars: individual (APM 08279+5255) 



1. Introduction 

Damped Lja systems (DLAs) comprise a class of the QSO absorption line systems 
characterized by high H I column density greater than 2 x 10^° [cm~^] (Wolfe et al. 1986). 
DLAs are one of the strongest probes to examine the evolution of metallicity and dust 
depletion in gaseous components of high redshift galaxies (e.g., Pettini 2004). Since there is 
a well-established correlation between A^Na i and A^hi (Ferlet et al. 1985; Bowen et al. 1995), 
Na I absorption, which has been observed extensively in the Galaxy, holds the potential to 
estimate A^h i for individual clouds in DLAs, which usually cannot be precisely measured 
because of heavy blending of H I absorption. Moreover, the column density ratio of Na I to 
Ca II provides a useful indicator of the degree of dust depletion of ISM clouds (e.g., Welsh et 
al. 1990; Crawford 1992; Bertin et al. 1993; Sembach et al. 1994; Price et al. 2001; Crawford 
et al. 2002) because sodium is hardly depleted from gas onto dust, while calcium can be 
heavily depleted (Savage & Sembach 1996). 

The main absorption lines of Na I are the optical AA 5891,5897 D doublet and the 
ultraviolet AA 3303.3, 3303.9 doublet. Observations exploiting the Na l D doublet hold the 
greater promise for comprehensive studies using Na I due to two important advantages: 



^ Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 
tory of Japan. 
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(1) the oscillator strength of the Na I D transitions are ~ 70 times greater than the Na I 
UV transitions (Morton 2003), and (2) Na I D absorption rarely overlaps with other metal 
absorption lines because the wavelengths of the Na I D doublet are significantly redder than 
most of other metal lines. Thus, Na I UV doublets can be detected only in the highest column 
density systems toward the brightest QSO. Weak Na I absorption systems are more likely to 
be detected with the Na I D doublet. However, observations of Na I D in high-redshift DLAs 
require near-infrared high resolution spectroscopy because the rest- frame optical wavelengths 
are redshifted into the near-infrared. The advent of high sensitivity and high resolution 
spectroscopy in the near-infrared with 8-meter class telescopes has enabled us to undertake 
such observations (e.g., Kobayashi et al. 2002, 2003, 2005). 

As a first trial for detecting Na I D absorption fines at high redshift {z > 1), we observed 
the extremely bright z^^ = 3.911 gravitationally-lensed QSO "APM 08279+5255" (Irwin et 
al. 1998). There are two dominant gravitationally-lensed images, "A" and "B," with a 
separation of 0'/38 and flux ratio of /b//a = 0.77 (Ibata et al. 1999). A third faint image, 
"C", was also reported by Ibata et al. between A and B, with a separation of A - C = 
0'.'15 and flux ratio fc/fA = 0.18. Pctitjean et al. (2000) reported probable faint Na I UV 
doublet absorption associated with the Zabs = 1-062 and 1.181 DLAs in optical HIRES (Vogt 
et al. 1994) spectra obtained by Ellison et al. (1999a,b). For both systems, Petitjean et al. 
concluded (1) log A'^hi ~ 21, (2) metalficity 0.3 — 1.0 Z0, (3) dust-to-metal ratio about half 
that of the Galaxy, and (4) a temperature of ~few hundred Kelvin. 

In this paper, we present near-infrared 1.18 -1.35 fim spectra of APM 08279+5255 with 
a resolution of about 50 km s^^. The spectra were obtained with the InfraRed Camera 
and Spectrograph (IRCS) using the adaptive optics (AO) system at the Japanese Subaru 
Telescope. The AO improved the FWHM of both images A and B such that the spectra of 
A and B were marginally separated. However, we combined the spectra of images A and 
B in order to maximize the signal-to-noise ratio and to compare our results with those by 
Petitjean et al. (2000). We will publish A-B separated spectra elsewhere. This paper is 
structured as follows: In §2, we describe the IRCS observations and data reduction of APM 
08279+5255. In §3, we present the Na I D absorption lines and the fitting results for the 
Zabs = 1-062, 1.173, 1.181 systems. We discuss our results and conclusions in §4. 



2. Observation &c Data Reduction 

Near-infrared 1.18-1.35 /im spectra of APM 08279+5255 were obtained on 2001 De- 
cember 25 UT, during the scientific commissioning run of the IRCS instrument. IRCS is 
a cassegrain mounted 0.9-5.5 /xm near-infrared camera and spectrograph (Tokunaga et al. 
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1998; Kobayashi et al. 2000) with an AO system (Takami et al. 2004) on the 8.2m Subaru 
Telescope (lye et al. 2004) located on Mauna Kea, Hawaii. 

The Subaru AO system uses a curvature wave front sensor with 36 control elements 
to compensate for a disturbed wavefront due to the earth's atmosphere. For an ~ 11 
wavefront reference star under excellent observing conditions, the AO provides a stellar 
image with a Strchl ratio of ~ 0.04 and a FWHM of ~ O'.'OT in the J-band (1.25 fj,m). The 

performance degrades with fainter reference stars. For these observations, we used APM 
08279+5255 itself {R ~ 15.6) as the wavefront reference star. We achieved a FWHM of 0'.'2 
in the J-band, which is much narrower than (0'.'5— 0'.'6) for typical seeing conditions at the 
Subaru Telescope. 

We used the cross-dispersed Echelle mode of the IRCS, which provides a spectral resolu- 
tion of 15 km s"^ (A/AA ~ 20000) with a 0'.'15 slit. The pixel scale is 0'.'075 (~ 7.5 km s'^) per 
pixel across the slit and 0'.'060 per pixel along the slilt. The entire J -band (1.18 — 1.35 /xm) 
was covered simultaneously in one IRCS exposure, with the cross-disperser in order 5 and 
the echelle in orders 42-48. Although the FWHM of the images were good throughout the 
observing time (about Qf!2-Qf.'5 in the J-band), we used a wider slit (0'/6) to maximize through- 
put. Thus, the spectral resolution was determined by the average image size ~ 0'/5, and the 
resulting spectral resolution was ~ 50 km s~^ (A/AA ~ 6000) in the observed wavelength 
range. 

The sht length was 3'/8, and the position angle was set to 32° so that both the lensed 
images A and B are in the slit. The telescope pointing was nodded along the slit by about 2" 
for each 600s exposure for sky-background and dark subtraction. Because of the fine pixel 
scale needed for AO images and the high spectral resolution, it is impossible to reach the 
background limit in 600s exposures, which are required to be less than the time variation of 
the sky OH emission. Eight sets of data were obtained, resulting in a total exposure time of 
9600s. Spectra of bright telluric standard stars (HD 92728: AOVs) at similar air mass were 
obtained in a similar fashion. 

All the data were reduced following standard procedures using the IRAF^ noao.imred. echelle 
package, including sky subtraction (subtraction of two frames), flat-fielding (halogen lamp 
with an integrating sphere), and aperture extraction. Argon lamp spectra that were taken at 
the end of the observing night were used for vacuum wavelength calibration. A Heliocentric 
correction of 9.64 km s~^ was applied to the vacuum wavelength calibrated spectra. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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3. Results 

In Figure 1, we present the 1.18 — 1.35 /xm spectrum of APM 08279+5255. The raw 
spectrum was smoothed with a 3 pixel (~ 20 km s~^) boxcar function for presentation. Na I 
D absorption lines from the Zabs = 1-062 and 1.181 DLAs are clearly seen. Several velocity 
components can be seen in the Na I D absorption lines from the Zabs = 1.181 DLA. Addition- 
ally Na I D absorption hnes associated with the Zabs = 1.173 Mg ll system (Petitjean et al. 
2000) are detected. This doublet absorption line was identified as Na I D lines because other 
metal hnes of other absorption systems hsted in Elhson et al. (1999a,b) cannot account for 
these two absorption features. 

We fit Voigt profiles to the Na I D absorption features assuming the component velocities 
identified by Petitjean et al. (2000) from metal absorption hnes in the optical spectra by 
Ellison et al. (1999a,b). The column density, Doppler width, and redshift of each component 
were evaluated using VPGUESS^ and VPFIT^ (Carswell et al. 1987). We fit two components 
for the Zabs = 1-062 system and five for the Zabs = 1-181 system. The absorption profiles 
are dominated by the instrument profile, which is almost identical to the Gaussian width 
FWHM ~ 50 km s'^ 

For the Zabs = 1-173 system, there is no a priori information on velocity components 
unlike the other systems- Also, there are several adjacent hot pixels on the right side of the 
Na I 5897 absorption line, which makes the fitting with multi velocity components difficult. 
Those hot pixels were difficult to remove because of their proximity to the Na I absorption 
line (see the feature labeled with x mark in Figure 4). Because of the above reasons, we 
performed only one component VPFIT (almost equivalent to a simple Gaussian fitting) for 
this system. The hot pixel regions were ignored for the fitting. In view of the statistical noise 
that can be estimated by point-to-point fiuctuations on the continuum, the single component 
fitting is almost satisfactory. 

The fitting results are summarized in Table 1. Figures 2,3,4 show the velocity profiles 
and fitting results for the Na I D hnes. Using atmospheric absorption lines in the object 

and standard spectra, we confirmed that the pixel shift along the dispersion direction was 
less than 1 pixel throughout the observing time. Therefore the systematic uncertainty of the 
velocity is less than 7.5 km s~^. Although the column densities and the Doppler widths could 



■^VPGUESS is a graphical interface to VPFIT written by Jochen Liske, 
http : / / www. cso . org /~] liskc / vpguess / 

''VPFIT is a Voigt profile fitting package provided by Robert F. Carswell, 
htt p : / / www . ast . cam. ac.uk/~rfc/vpfit. html 
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not be directly determined because of insufficient instrumental spectral resolution, they can 
also be constrained by the doublet ratio [i.e., equivalent width ratio W{5891)/W{5S97)]. 
As the doublet ratio approaches unity, uncertainty in column density increases because the 
column density and Doppler width become degenerate for a fixed equivalent width. 

Because the physical distance (~ 1 kpc) between A and B at 2; ~ 1 is much larger than 
the typical spatial scale of Na I clouds in our Galaxy (a few tens to several thousands AU e.g.. 
Price et al. 2001), the images A and B are most likely covered unequally. This introduces an 
additional systematic uncertainty in the inferred properties of the Na I clouds. An estimation 
of the systematic uncertainty depends on the saturation level of the Na I absorption lines 
in the A+B combined spectrum. It can range from ~0.3 dex for a moderately saturated 
case to ~ 1 dex for a heavily saturated case. The estimated systematic uncertainty for each 
component is shown in Table 1. 

Because the resolution of the optical spectrum is higher than our near-infrared spec- 
trum, the estimated redshift of each component by Petitjean et al. (2000) should be more 
accurate than ours. Hereafter, we use the velocities reported by Petitjean et al. (2000) 
for our discussion of each component, except for the newly identified components from our 
spectrum. 

3.1. DLA at Zabs = 1-062 

In this system, we found two Na I components at -|-53 and -1-108 km s~^, which we 
identify as the components at -|-55 and -1-120 km s"-*^ by Petitjean et al. (2000). While 
Petitjean et al. found Na I component only at the -|-55 km s~^, we newly identified Na I 
component at the -1-120 km s~^. The estimated column densities are log Nnsi = 12.9 and 11.7, 
respectively. These are consistent with Petitjean et al., who detected the former component 
with the similar column density, but not the latter component to an upper- limit of log ^Nai < 
12.8. 



3.2. DLA at Zabs = 1.181 

In this system, we found five Na I components. We identified the three components at 
—99, —75, and +48 km s~^ as the components at —105, —80, and -|-40 km s~^ reported by 
Petitjean et al. (2000). While Petitjean et al. found Na I component only at the —80 km s~^, 
we newly identified Na I components at the —105, and +40 km s~^. We identified further two 
new Na I components at +4, and +101 km s~^, which were not identified in Petitjean et al 
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probably due to contamination by other absorption lines. The estimated column density of 
the —80 km s~^ component is log NNai = 12.5, which is not consistent with Petitjean et al. 
However, our evaluation of the column density of —80 km s~^ component is not so robust 
because Mg I A2852 absorption associated with an Mg II absorption system at Zabs = 3.502 
(Kondo et al. in preparation) overlaps with the Na I D A5891. The estimated column density 
of the —105 km s~^ component is log NNai = 12.0, which is consistent with the upper-limit 
of log NNai < 12.4 found by Petitjean et al. 

3.3. Mg II system at Zabs = 1-173 

In this system, we found a Na I component at +204 km s^"*^. Petitjean et al. (2000) 
did not find any metal absorption lines at this velocity. The detected Na I absorption line 
probably arises in a cloud that was previously not in front of the quasar at this velocity. 
We discuss the detail of this newly identified component in §4.2. We could not detect Na I 
absorption lines associated with the four components around km s~^ reported by Petitjean 
et al., probably because of their low Na I column density. 

4. Discussion 

4.1. Hi Column Densities of Each Velocity Component 

Following Petitjean et al. (2000), we estimated Nhi for all components using the em- 
pirical relationship, log Nhi = 0.688 log Nwai + 12.16 (11.6 < logNwai < 13.4) (Sembach et 
al. 1993; Diplas & Savage 1994; Bowen et al. 1995) established for the Galactic interstellar 
absorption clouds. In Table 2, we summarize our results for Nhi- While this empirical rela- 
tionship does not hold for log NNai < H (Welty et al. 1994; Wakker & Mathis 2000; Wakker 
2001) nor for log NNai > 12.9 (Ferlet et al. 1985), the column densities of all the detected 
components are within the appropriate range of this relationship. The estimated Nhi in 
various components, however, is subject to systematic uncertainties related to the following 
issues: 

1. Low metallicity 

It is likely that the column density ratio NNai/Nni is in proportion to the metallicity. 
Vladilo et al. (1993) found that the column density ratio in the Large Magellanic Cloud 
(LMC) decreases in proportion to the metallicity (~ 0.3 Zq, Peimbert et al. 1974; Pagel 
et al. 1977). Because the metallicities of the Zabs = 1.062 and 1.181 DLAs are estimated 
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to be in the range of 0.3 — 1.0 Zq (Petitjean et al. 2000), it is appropriate to interpret 
the estimated Nhi in these systems as lower-hmits. 

2. Contamination from a metal absorption hne of other redshift 

Wc could not evaluate the reliable column density of the —80 km component in 
the Zabs = 1.181 system because Mg I A2852 absorption associated with a strong 
Mg II system at Zabs = 3.502 (Kondo et al. in preparation) is blended with the Na I D 
A 5891 absorption. Assuming the correlation, log NMgi/NMgii = —0.73 log NMgii + 7.6 
(Churchill et al. 2003b) and the same Doppler width (b = 13.7 km s~^) as the Mg ll 
absorption lines, we estimated the equivalent width of the Mg I A2852 absorption. The 
H I column density of this component was estimated after removing this contamina- 
tion from Mg I, and this introduced systematic uncertainty. Even after removing the 
contamination, there still is another large uncertainty, because the resultant doublet 
ratio of the Na I absorption lines is close to unity. 

3. Coverage of gravitationally-lensed images A and B 

As stated above, it is likely that the Na l clouds do not equally cover the gravitationally- 
lensed images A and B. Thus, the estimated Na I column density is uncertain as 
described in §3. The estimated Nhi varies accordingly. 

While Petitjean et al. (2000) suggested that there are 40 and 20 small clouds in the 
Zabs = 1-062, and 1.181 systems, respectively, our results show that there are at least 2 and 
5 high column density clouds in each system. Most interestingly, we find that each cloud in 
the Zabs = 1.181 system has a H I column density that is classified as a DLA. This is the 
first indication of an absorber in which the individual clouds are each of such high column 
density. This is in contrast to the findings of Churchill et al. (2003a), who constrained the 
size of DLA clouds be less than 25 pc in the lensed QSO 0957+56lA,B and inferred that 
only one of the clouds gave rise to the DLA H I column density. This would suggest that 
the Zabs = 1.181 DLA, having log A^hi ^ 21.2, may be probing a very massive and large 
underlying structure (see, for example, Turnshek et al. 2004). 

4.2. Detection of a New DLAs at Zabs = 1173: Proper Motion of a cloud of a 

few 100 AU? 

We detected a new component at +204 km s~^ associated with the Zabs = 1.173 Mg ll 
absorption system. Based upon scaling from the high Na I column density, this system is 
probably a DLA rather than a Lyman-limit Mg ll system (Petitjean et al. 2000). Therefore, 
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this system might contribute to the gravitational- lensing of the QSO, together with the other 
three DLAs, at Zabs = 1-062, 1.181, and 2.974, as suggested by Petitjean et al. Despite the 
large Na I column density, Petitjean et al. did not identify any metal lines for this velocity 
component in the optical Keck spectrum. Similarly, in a subsequent search, we also could 
not identify any Mg II, Mg I, and Ca ll absorption lines for this component in the HIRES 
spectrum^ (see Figure 4). 

Because there is a time-span of 3.6 years between our observation and the optical obser- 
vation in the local-frame, or ~1.7 years in the rest-frame (zabs = 1-173) considering the effect 
of cosmological time dilation by a factor of (1 + Zabs)~^, it is probable that the newly identi- 
fied cloud appeared in the line of sight toward the QSO within that time-span. If the cloud 
is moving in a tangential direction in the sky with a velocity of ~ 400 km s~^, the travel 
distance is ~ 200 AU in ~1.7 years in the rest- frame. The assumed velocity is comparable to 
the combined velocity of the motion of cloud in high-redshift galaxies (~ 50 to ~ 260 km s~^, 
Erb et al. 2003) and the proper motion of field galaxies which should be smaller than the 
typical velocity dispersion of the cluster of galaxies (e.g., between ~ 500 and 1000 km s~^. 
Gal & Lubin 2004). This sets the upper-limit of the cloud size and it is consistent with the 
size of the Na I clouds in the Galaxy, which ranges from a few to several thousands AU 
(e.g., Meyer et al. 1996; Meyer & Lauroesch 1999; Welty & Fitzpatrick 2001; Andrews et al. 
2001; Lauroesch et. al. 2000, 2003). Our result shows that we can catch the very small-scale 
structures of the cold (~ 100 K) neutral gas component in high-z QSO absorption systems 
with time variation of Na I absorption features. The observable spatial scale with this tech- 
nique is much smaller than that with the observation of gravitationally-lensed QSOs (see 
discussion in Ellison et al. 2004). 



4.3. Dust Depletion and Chemical Uniformity in DLAs 

The column density ratio NNai/Ncaii has often been used as a good indicator of the 
degree of dust depletion because sodium is hardly depleted from gas onto dust grains while 
calcium can be heavily depleted (Savage & Sembach 1996). In Table 3, we hst the ratio 

NNai/Ncaii of cach component. Ncaii for the -|-120 km s~^ component at Zabs = 1-062 and 
for the —105 and —80 km s~^ components at Zabs = 1-181 arc taken from Petitjean et al. 
(2000), while those for the +4, +40, and -1-101 km s~^ components at Zabs = 1-181 were 



5ftp://ftp-ast-cam.ac.uk/pub/papers/APM08279 



-10- 



evaluated^ using archived optical spectra (Ellison et al. 1999a,b). 

We estimated the 3a upper-hmit of Ncaii for the +4, and +101 km s~^ components at 
Zabs = 1.181, because no Ca ll absorption line was detected in the optical spectra. In Figure 
4, we show a log [NNai/Ncaii] - log NnoI plot for the DLAs as well as Na l clouds in the 
Galaxy (Vallerga et al. 1993) and in the LMC (Vladilo et al. 1993). 

We found that log [NNai/Ncaii] in these DLAs is systematically smaller than those in the 
Galaxy but similar to those in the LMC (Welty et al. 1999), suggesting that dust depletion 
is generally smaller in a lower metallicity environment. Our results are consistent with the 
results of Vladilo (1998), who shown that dust-to-metal ratios of DLAs are similar to that 
of the LMC but smaller than that of the Galaxy, using the column density ratio of zinc to 
chromium. 

It is also interesting to consider the uniformity in the NNai/Ncair ratio across the ~ 200 
km s~^ velocity spread of the Zabs = 1.181 DLA. As presented in Table 3, it would appear that 
the ratios are consistent with being uniform with component velocity, though there is large 
uncertainty in this statement. If so, this is consistent with the findings of Prochaska (2003a) 
and Churchill et al. (2003a), who reported uniform column density ratios as a function of 
velocity in DLAs. It is also found by Churchill et al. that the uniformity is spatial, at least 
on the scale of 200 pc. If the NNai/Ncaii ratio is more-or-less constant across the profiles 
in the Zabs = 1.181 DLA, it would suggest that the chemical, ionization, and dust depletion 
levels are uniform across the clouds. 



4.4. The Zabs = 1181 System: Cold Dense Gas Clouds near the Center of a 

Galeixy ? 

All five components of the Zabs = 1.181 system have log [NNai/Ncaii] > 0, which is a 
characteristic of cold dense gas in the Galaxy (known as the "the Routly-Spitzer effect", 
e.g., Routly & Spitzer 1952; Vallerga et al. 1993; Sembach et al. 1994). Because there are 
five cold dense and very high column density clouds for the Zabs = 1.181 system, we suggest 
that the line of sight may be passing near the center of the host galaxy. The relatively high 
metallicity of this system (0.3 — 1.0 Zq, Petitjean et al. 2000) compared to that typical of z~l 
DLAs (~ 0.1 Zq, Prochaska et al. 2003b) also supports the above idea, because the metallicity 



^We evaluated Ca ii column density for those systems using VPFIT. We evaluated Ca ii column den- 
sities also for the +120 km s~^ component at Zabs = 1.062 and the —105, and —80 km s~^ components at 
Zabs = 1.181 and confirmed that our fitting results are consistent with those estimated by Petitjean et al. 
(2000). 
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of DLA absorption is higher in the inner region galaxies because of a metaUicity gradient 
(Chen et al. 2005). Therefore, the host galaxy of this system might be the most significant 
contributor to the gravitational-lens among four DLAs toward APM 08279+5255. Note that 
some DLAs were observed to have an impact parameter larger than 65 kpc (Churchill et al. 
2005) and above argument should be taken with caution. 

We are grateful to all of the IRCS and AO team members and the Subaru Telescope 
observing staffs for their efforts, which made it possible to obtain these data. Wc arc also 
grateful to Ellison ct al. for making their excellent data available in the Internet and to 
Jochen Liske for kindly answering our questions about VPGUESS. We thank the National 
Astronomical Observatory of Japan for their financial support and encouragement for the 
construction of the IRCS and AO. Y.M. is financially supported by the Japan Society for 
the Promotion of Science (JSPS). 
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Fig. 1. — Detected Na I D absorption lines at Zabs = 1-062 (left panel) and Zabs = 1-173, 1.181 
(right panel). Vertical axis is normalized flux to the continuum level and horizontal axis 
is local-frame wavelength in both panels. The spectrum was smoothed with a 3 pixel 
(~ 20 km s~^) boxcar function. Longer hues indicate Na I D A5891, while shorter hues 
indicate Na I D A5897. The left spectrum is noisier than the right one because the Na I D 
absorption lines at Zabs = 1-062 were located on the edge of the detector where the detec- 
tor noise is significantly larger. The bottom panel shows the telluric absorption spectrum 
calculated with ATRAN package software (Lord 1992). 
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Fig. 2. — Velocity profile of metal lines in Zabs = 1.062 system. Horizontal axis shows the 
heliocentric velocity from z=l. 06230. Vertical axis is normalized flux to the continuum level. 
Top two panels show the velocity profiles of Na I D AA5891, 5897 absorption lines obtained 
with IRCS with the velocity resolution of ~ 50 km s~^. Bottom four panels show the velocity 
profile of Na I A3303.9, Ca ll A3934, Mg ll, A2803 and Ti ll A3242 absorption lines obtained 
with HIRES with the velocity resolution of ~ 8 km s~^. Vertical dashed lines mark the 
velocity components found by Petitjean et al. (2000). The blue line shows the Na I velocity 
component which has been detected with Na I AA3303.3, 3303.9 (Petitjean et al. 2000) and 
the red line shows the newly detected Na I velocity component determined from the VPFIT 
program, which was also seen in the Ca ll, Mg ll, and Ti ll line profiles. Gray lines show the 
sum of all components. In the HIRES spectra, unmarked absorption lines are either from 
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Fig. 3. — Velocity profile of metal lines in the z^bs = 1-18 system. Notations are same 
as those in Figure 2. The horizontal axis shows the heliocentric velocity from z= 1.1807. 
Bottom four panels show the velocity profile of Na I A3303.9, Ca ll A3934, Mg ll A2803 
and Mg I A2825 absorption lines obtained with HIRES. Vertical dashed lines mark veloc- 
ity components (V = —105, —80, +40 km s~^) found by Petitjean et al. (2000) and those 
(V = +4, +101 km s~-^) newly found with Na I D data. 



- 18 - 



1 

0.95 
0.9 

1 

0.95 
0.9 

1 

0.9 

6 0.7 
U 

o 

IZi 



h NalD 
X5891 



I — I 

0) 
N 



- NalD 
X5897 



1 - 



0.5 


1 



0.5 - 



-200 



Call 
X3934 





X 



X 



200 



400 



Heliocentric Velocity (km s ^) 



Fig. 4. — Velocity profile of metal lines form Zabs = 1.173 system. Notation are same 
as those in Figure 2. Horizontal axis shows the heliocentric velocity from z=1.1727. 
Bottom three panels show the velocity profiles of Ca ll A3934, Mg ll A2803 and 
Mg I A2852 absorption lines obtained with HIRES. Vertical dashed lines mark veloc- 
ity components (V = — 20, 0, +35, +50 km s~^) found by Petitjean et al. (2000) and those 
(V = +204 km s~^) newly found with Na I D data. Spurious features labeled with x mark 
on the right side of A5897 absorption are due to hot pixels. See the main text for detail. 
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Fig. 5. — The relation of log NNai and log [NNai/Ncaii]- Red points show our DLA data for 
Zabs = 1.181. A red circle shows our DLA data for Zabs = 1.062. Black and blue points 
show the data for our Galaxy (Vallerga et al. 1993) and the LMC (Vladilo et al. 1993) , 
respectively. The uncertainty of those data is less than 20%. Red filled square shows the 
data for the V = —80 km s~^ velocity component at Zabs = 1-181, whose uncertainty is not 
shown because it could not be accurately estimated due to the accidental overlap of the 
Zabs = 3.502 Mg I absorption line (see main text). 
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Table 1. The velocity components of the detected Na I D absorption hnes 



Our results Petitjean et al. (2000) 



system 


Z^abs 


Ya,b 


log NNal 


b 


ya.c 


logNNal 


b 






(km s""'^) 




(km s~^) 


(km s-i) 




(km s^^) 


1.062 


1.06266±0.00001 


+53±2 


12. g'^ ±1.6 


1.1±0.7 


+55 


12.9 


1.5 


1.062 


1.06305±0.00002 


+108±3 


11. 7^= ±0.1 


7±12 


+120 


< 12.8 


1.5 


1.173 


1.17418±0.00001 


+204±1 


12.3<* ±0.6 


1.0±0.4 








1.181 


1.17998±0.00006 


-99±8 


12.0^1 ±0.5 


1.3±2 


-105 


< 12.4 


2.5 


1.181 


1.18015±0.00003 


-75±4 


12.5'' ±0.1 


4±2 


-80 


13.5 


0.8 


1.181 


1.18073±0.00005 


+4±6 


11.7' ±0.5 


0.8±1 








1.181 


1.18105±0.00003 


+48±4 


12.0** ±0.1 


3.6±3 


+40 






1.181 


1.18143±0.00002 


+101±2 


12.2«^ ±0.2 


1.8±0.6 









'^The velocity zero point is set to z=l. 06230, 1.17270, 1.18070 for each system. 

''There is a systematic uncertainty which is less than ~ 7.5 km s~^. See main text for detail. 

■^We estimated the velocity with an accuracy of 5 km s"'^ from Figure 7,8 of Petitjean et al. (2000) by eyes 
because Petitjean et aL (2000) do not Ust the velocities. 

''There is a systematic uncertainty of about 1 dex because the gravitationally-lensed images are not 
resolved. See main text for detail. 

''There is an additional systematic uncertainty of about 0.3 dex because the gravitationally-lensed images 
are not resolved. See main text for detail. 

There is a possible significant systematic uncertainty because of the accidental overlap of the Zg^g = 3.502 
Mg I absorption line. See main text for detail. 
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Table 2. Estimated Nhi for the Na I D absorption lines 







Our 


results 


Petitjean et al. (2000) 


system 


V 


log NNal 


log Nhi 


log NNal 


log Nhi 


1.062 




12.9"^ ±1.6 


21.0^ ±1.1 


12.91 ±0.04 


21.0 ±0.0 


1.062 


+ 120^ 


11.6" ±0.1 


20.2" ±0.1 


< 12.8 


< 21.0 


1.173 


+204 


12.3'^ ±0.6 


20.6^' ±0.4 






1.181 


-105^ 


12.0'^ ±0.5 


20.4^ ±0.3 


< 12.40 


< 20.7 


1.181 


-80^ 


12.9^^ 


20.7'^ 


13.5 ±0.11 


21.4 ±0.1 


1.181 


+4 


11.7" ±0.5 


20.2" ±0.3 






1.181 


+40^ 


12.0" ±0.1 


20.4" ±0.1 






1.181 


+101 


12.2" ±0.2 


20.5" ±0.1 







^Velocities from figures in Petitjean et al. (2000). See note b in Table 
1 for detail. 

''There is an additional systematic uncertainty of about 1 dex. See 
Table 1. 

"There is an additional systematic uncertainty of about 0.3 dex. See 
Table 1. 

'^There is a possible significant systematic uncertainty. The value 
shown here is after removing of the contamination from the Mg I ab- 
sorption at Zabs = 3.502. 
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Table 3. Comparison of log NnsI and log [NNai/Ncaii] for the Na I D absorption lines 



syst cm 


V 


log NNal 


log Ncaii 


log[XKal/N(.a,Il] 


1.062 


+50 


12.9^+1.6 


. . A 


. . .d 


1.062 


+110 


11.7'^ ±0.1 


11.73 ± 0.03 


-0.02 ± 0.1 


1.181 


-100 


12.0^ ±0.5 


11.33 ±0.39 


0.7 ±0.6 


1.181 


-75 


12. 9^= 


11.79 ±0.35 


1.1 


1.181 


+4 


11.7^ ±0.5 


< 11 


> 0.7 


1.181 


+40 


12. O*' ±0.1 


11.4 ±0.2 


0.6 ±0.1 


1.181 


+101 


12.2'" ±0.2 


< 11 


> 1.2 



^There is an additional systematic uncertainty of about 1 
dex. See Table 1. 

^There is an additional systematic uncertainty of about 0.3 
dex. See Table 1. 

'^There is a possibly significant systematic uncertainty. The 
value shown here is after removing the contamination from the 
Mg I absorption at Zabs = 

3.502. 

^No measurement was possible due to the blending effects. 
See the note of TABLE 2 in Petitjean et al. (2000). 



